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Abstract. We measure the two-point redshift -space cor- 
relation function of loose groups of galaxies, ^gg(s) , for 
the combined CfA2 and SSRS2 surveys. Our combined 
group catalog constitutes the largest homogeneous sam- 
ple available (885 groups). We compare ^gg(s) with the 
correlation function of galaxies, ^gg(s) , in the same vol- 
ume. We find that groups are significantly more clustered 
than galaxies: < <^GG/Cgg > =1.64 ±0.16. A similar result 
holds when we analyze a volume-limited sample (distance 
limit 78 h~^ Mpc) of 139 groups. For these groups, with 
median velocity dispersion ^ 200 km s~^ and mean 
group separation d ■^^ 16 h~^ Mpc, we find that the corre- 
lation length is So = 8 ± 1 h~^ Mpc, which is significantly 
smaller than that found for rich clusters. We conclude that 
clustering properties of loose groups of galaxies are inter- 
mediate between galaxies and rich clusters. Moreover, we 
find evidence that group clustering depends on physical 
properties of groups: correlation strengthens for increas- 
ing (Tv 

Key vi^ords: cosmology: large-scale structure of Universe 
- galaxies: clusters: general - galaxies: statistics 



1. Introduction 

Loose groups of galaxies, the low-mass tail of the mass dis- 
tribution of galaxy systems, fill an important gap in the 
mass range from galaxies to rich clusters. Until now, clus- 
tering properties of loose groups have been studied on the 
basis of rather small samples. The results are consequently 
uncertain and even contradictory. Nevertheless, cluster- 
ing properties of groups are shown to be robust against 
the choice of the identification algorithm, provided sys- 
tems are identified with comparable number overdensity 



thresholds (Frederic 1995b). 

The two-point correlation function, CF, of galaxies 
and of galaxy systems constitutes an important measure 



of the l arge-s cale distribution of galaxies (e.g 
Peebles |1983[ Bahcall 
Tucker et al 



1988: 



Davis & 

Soneira 1983; de Lapparent et 
Croft et al. |I997D. From 



1997 



al 

galaxies to clusters the two-point correlation function in 
redshift-space, ^(s), (e.g., Peebles 1980), is consistent, 
within errors, with a power-law form ^(s) = (s/sq)"'^ 
with 7 ~ 1.5 — 2 for a variety of systems. The correla- 
tion length. So, ranges from about 5-7.5 h^^ Mpc for 
galaxies (e.g., Davis & Peebles 1983|; L oveday et al. |1996 
Tucker et al. 
to So 15 h 



1997| ; Willmer et al. |1998| ; Guzzo et al. [l999^ 
^ Mpc for galaxy clusters (e.g., Bahcall & 



Soneira 1983; Postman et al. 1991; Peacock & West 



Miller et al. 1999 



1992 



Croft et al. |1997|; Abadi et al. [1998^ Bor gani et al 



199£) 



1999 



Moscardini et al 
As far as loose groups are concerned, previous deter- 
minations of the CF are very uncertain. From the study of 
137 groups (wit hin Cf Al) and 87 groups (withi n SSR Sl), 
Jing & Zhang ( |l988| ) and Maia & da Costa (|l990D re- 
spectively find that the group-group CF, ^gg(s) , has a 
lower amplitude than the galaxy-galaxy CF, Cgg(s) ■ An- 
alyzing 128 groups in a sub-volume of CfA2N, Ramella 
et al. (1990; hereafter RGH90) find that the amplitudes 
of ^gg{s ) and Cgg(s) are consistent (see also Kalin kov & : 
Kuneva |1990D . Finally, Trasarti-Battistoni et al. ( |1997D 
study 192 groups in the Perseus-Pisces region and find 
that the amplitude of ^gg{s) exceeds that of ^gg{s) . 

The theoretical expectations for the relative strength 
of the group and galaxy clustering are also contradictory. 
Frederic (1995b) determines the correlation function for 
galaxy and group halos in CDM numerical simulations by 



Gelb (1992) and finds that groups are more strongly cor- 
related than galaxies. In contrast, Kashlinsky (1987), on 
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the basis of an analytical approach to the clustering prop- 
erties of collapsed systems of different masses, concludes 
that groups and individual galaxies should be correlated 
with the same amplitude. 

Here we compute the two-point correlation function 
(in redshift space) for 885 groups of galaxies identified 
in the combined CfA2 and SSRS2 redshift surveys. This 
sample is characterized by its large extent (more than five 
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Table 1. Data Samples 



Sample 


Ng 




CfA2N 


395 


5426 


CfA2S 


139 


2104 


SSRS2N 


123 


1697 


SSRS2S 


228 


3083 


TOTAL 


885 


12290 



times the volumes previously studied) and by the homo- 
geneity of the identification process (the friends-of~friends 
algorithm FOFA; Ramella ct al. [1997| - hereafter RPG97). 
Moreover, we compare the group-group CF to that com- 
puted for galaxies in order to determine the relative clus- 
tering properties of groups and galaxies. Because we use 
the same galaxy sample where groups are identified, we 
avoid possible effects of fluctuations due to the volume 
sampled. 

In Sect. 2 we briefly describe the data; in Sect. 3 we 
describe the estimation of the two-point correlation func- 
tion; in Sect. 4 we compute the correlation function of 
groups and compare it to that for galaxies; in Sect. 5 we 
summarize our results and draw our conclusions. 

Throughout the paper, errors are at the 68% confi- 
dence level, and the Hubble constant is Hq = 100 h Mpc^^ 
km s"""^ . 

2. Galaxy and groups catalogs 

We extract the sample of galaxies from the CfA2 North 



(CfA2N) and South (CfA2S) ( Geller & Huchra |1989; 
Huchra et al. |1995| ; Falco et al. |1999D , and the SSRS2 
North (SSRS2N) and South (SSRS2S) (da Costa et al. 
1998| ) redshift surveys. These surveys are complete to 
mB(o) — 15.5 and cover more than one-third of the sky, 
i.e. most of the extragalactic sky. The original papers con- 
tain detailed descriptions of the observations and of the 
data reduction. The velocities we use are heliocentric; they 
include corrections for solar motions with respect to the 
Local Group and for infall toward the center of the Virgo 
cluster (see RPG97 for details). As in previous analyses 
of th e CfA2 surveys (e.g., Park et al. [1994 Marzke et al. 
1995| ), we discard regions of large galactic extinction. The 



total sample includes 13435 galaxies with radial velocity 
V < 15000 km s-i . 

We use the catalogs of groups identified within CfA2N 
by RPG97 and within SSRS2 by Ramella et al. (in prepa- 
ration). The identification method is a friends-of-friends 
algorithm (FOFA) which selects systems of at least three 
members above a given number density threshold in red- 
shift space. In particular, RPG97 and Ramella et al. (in 
preparation) use the number density threshold <5pn/pn = 
80 and a line-of-sight link Vq = 350 km s~^ at the fidu- 
cial velocity V{ = 1000 km s~^ . We run FOFA with these 
parameters on CfA2S and produce a group catalog for 
this survey, too. The combined catalog contains a total of 
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Fig. 1. We show the projection on the sky of the group 
and galaxy samples (upper and lower panel, respectively). 

885 groups that constitute a homogeneous set of systems 
objectively identified in redshift space. 

The group catalogs are limited to radial velocities V < 
12000 km s^^ , but members are allowed out toV< 15000 
km s~^ . We confine the galaxy sample to V < 12000 
km s~^ and are left with a total of 12290 galaxies. 

Table ^ lists the numbers of groups, Nq, and the num- 
bers of galaxies, iVg, for each sample. Fig. |l| shows the 
distribution of the galaxy and group samples on the sky. 

3. Estimation of the correlation function 

We compute the two-point correlation functions in red- 
shift space for groups and galaxies (hereafter ^gg(s) , and 
^gg(s) , respectively). The formalism in the two cases is 
the same. We define the separation in the redshift space, 
s, as : 



Ho 



(1) 



where Vi and Vj are the velocities of two groups (or galax- 
ies) separa ted by an angle 9ij on the sky. Following Hamil- 
ton (1993) we estimate ^(s) with: 



^ DD{s)RR{s) 



(2) 



where DD{s), RR{s), and DR{s) are the number of data- 
data, random-random, and data-random pairs, with sep- 
arations in the interval (s, s + ds). We build the control 
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sample by filling the survey volume with a uniform ran- 
dom distribution of the same number of points as in the 
data. The points are distributed in depth according to the 
selection function of the surveys, ^{V). 

In order to decrease the statistical fluctuations in the 
determination of ^(s), we average the results obtained us- 
ing several different realizations of the control sample. We 
compute 50 realizations in the case of groups and 5 in the 
case of galaxies. 

Unless otherwise specified, we compute the "weighted" 
correlation function by substituting the counts of pairs 



with '. 



the weighted sum of pairs, which takes into 



account the selection effects of the sample used. In the case 
of a sample characterized by the same selection function, 
volumes are equally weighted and Wi = l/<i>(T4) is 
the weight of a group (or galaxy) with velocity Vi. The 
appropriate selection function fo r a m agnitude-li mited 
sample (e.g., de La pparent et al 
Willmer et al. |1998| ) is: 



1988; Park et al. 1994 



J— CX 



(l){M)dM 



(3) 



where 4'{M) is the Schechter (1976) form of the lumi- 
nosity function. Afmax is a low luminosity cut-off. We 
chose Minax — —14.5, the absolute magnitude correspond- 
ing to the limiting apparent magnitude of the survey at 
the fiducial velocity Vf — 1000 km s^^ . This value of 
Vf is the same as in RPG97. The Schechter parameters 
of the galaxy luminosity function, before the Malmquist 
bias correction, are: M* = —19.1, a = —1.1 for CfA2, and 
M* = -19.7, a = -1.2 for SSRS2 (Marzke et al. |l994; 



1998). 



We assume that the group selection function is the 
same as for galaxies. In fact, the velocity distributions 
of groups, A'^g(^^), and of galaxies, Ng{V), are not sig- 
nificantly different according to the Kolmogorov-Smirnov 
t est ( d. Fig. I). RG II90 an d Trasarti-Battistoni et al. 
( 119971) , and Frederic ( |l995b| ) make the same assumption 
for observed and simulated catalogs, respectively. 

The different luminosity functions of CfA2 and SSRS2 
correspond to different selection functions. For this reason 
we assign to a group (galaxy) i, belonging to the subsam- 
ple fc, the weight given by: 



1 



Wi 



*k(Vi)nk' 



(4) 



where rik is the mean number density of groups (galaxies) 
of that subsample (e.g. Hermit et al. |1996 ). We compute 
the density as rik = 1/Vk • Si[l/<I'(Vi)], where the sum is 
over all the groups (galaxies) of the subsample volume. 



Vk (Yahil et al. 1991). In our analysis, Vk is the effective 



volume of the subsample. Because the different selection 
functions, we also build a control sample for each sub- 
sample separately, and, conservatively, we do not consider 
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Fig. 2. The velocity distributions of groups (solid line) 
and of galaxies (dotted line) are compared. The error 
bands represent 1 — a Poissonian errors. 




0.1 



Fig. 3. The group-group (filled circles) , and galaxy- 
galaxy (open circles) CFs, as computed for the combined 
sample. The error bands represent 1 — cr bootstrap errors. 



pairs of groups (galaxies) linking two different subsam- 
ples. In this way we also avoid crossing large unsurveyed 
regions of the sky. 

We compute the errors on ^(s) from 100 bootstrap re- 
samplings of the data (e.g.. Mo et al. 1992| ). Note that the 
bootstrap-resampling technique, which overestimates the 
error in individual bins, represents a conservative choice 
in this work. 



4. The group— group correlation function 

We plot the group-group CF, ^gg(s) , in Fig. |3[ In the 
same figure we also plot the galaxy-galaxy CF, Cgg(s) . 
On small scales (s <^ 3.5 Mpc ) ^gg(s) starts dropping 
because of the anti-correlation due to the typical size of 
groups. On large scales (s ^ 15 h^^ Mpc) the signal-to- 
noise ratio of ^00(5) drops drastically. We thus limit our 
analysis to the separation range 3.5 ^ s ^ 15 h~^ Mpc. 
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Fig. 4. The ratio ^gg /Cgg j where the CFs are those com- 
puted for the combined sample in Fig. 3. 




Fig. 5. The ratio ^cc/Cgg fo^' the combined sample (thin 
line), but, now, ^gg{s) is computed also for groups with 
at least five member galaxies (thick line) and for groups 
with less than five member galaxies (dashed line). 



The main physical result in Fig. || is that ^gg (s) has 
a larger amplitude than £,gg{s) . This property of the CFs 
is also evident in Fig. ^ where we plot the ratio £,gg{s) / 
Cgg(s) on a linear scale. Over the s-range of interest, the 
values of the ratio are roughly constant within the errors. 
In order to give an estimate of the relative behavior of 
groups and galaxies we compute the mean of the values of 
the ratio. We obtain < ^Gc/^g > =1.64 ± 0.16. 



4-.1. The CF of rich groups 

Groups with a number of members A'mcm > 5 are gener- 
ally reliable, as shown both by optical and X-ray analyses 



(Ramella et al. 1995; Mahdavi et al. 1997). On the other 



hand, the reliability of groups with fewer members is of- 
ten questionable. In particular, the analysis of the CfA2N 
survey performed by RPG 97, and the analysis of a CDM 
model by Frederic (1995a) show that a significant fraction 



of the triples and quadruples in group catalogs could be 
spurious. 

We consider the 321 rich groups with Nn-icm > 5, and 
find again that groups are more correlated than galax- 
ies. Moreover, we find evidence that the CF computed 
for richer groups is higher than the CF computed for 
poorer groups, i.e. < CGG,rich/^gg >= 1-48 ± 0.16 and 

< eGG,poor/egg >= 1.10 ± 0.14 (Cf. Fig. I). 

The component of spurious groups among poor groups 
could be responsible for the lower amplitude of the 
fGG(s) of poor groups compared to the £,GGis) of rich 
groups. In fact, spurious groups should be distributed like 
non-member galaxies. However, at least part of the ob- 
served higher clustering amplitude of rich groups could 
be due to the existence of a clustering amplitude vs rich- 
ness relationship. The relationship has been discussed for 
a variety of systems by several authors (e.g. Bahcall 
Croft et al. 



West 1992 



1997; Miller et al. 1999). Richness 



is usually taken as a measure of the mass of the system, 
mass being the real, interesting physical quantity directly 
related to the predictions of cosmological models. 

4-2. The CF in the volume-limited sample 

For our groups, richness is not a good physical parame- 
ter. A better parameter is the group (line-of-sight) veloc- 
ity dispersion, ctv (e.g. RPG97). In a magnitude-limited 
sample any group selection based on velocity dispersion 
will affect the selection function in an "a priori" un- 
known way. To avoid this problem, we analyze the volume- 
limited group sample built by Ramella et al. (in prepa- 
ration) who run an appropriately modified version of 
FOFA within volume-limited sub-samples of the CfA2 
and SSRS2 galaxy surveys. In particular, we consider the 
139 distance limited groups within V < 7800 km s~^ , 
roughly corresponding to the effective depth of CfA2. We 
cut the volume-limited galaxy catalogs in the same way. 
Within this sample we compute the "unweighted" CF es- 
timator, i.e. we set u; = 1 for all groups/galaxies. For this 
sample the useful s-range is 3.5 ^ s < 12 h^^ Mpc (see 
Fig. I). 

We find that the ratio < ^GG/'?gg > of the total 
volume-limited sample is < ^GG/Cgg > = 1.58±0.10, sim- 
ilar to that computed for the magnitude-limited sample 
(< ^GG/Cgg > ^ 1-6). This result reassures us about the 
reliability of the selection function we assume for groups. 

In order to check a possible dependence of ^gg(s) on 
(Tv, we divide the group volume-limited sample into two 
subsamples of equal size, one subsample containing groups 
with Ctv > 214 km s~^ , the other including the remain- 
ing low velocity dispersion groups. We find that high- 
CTv systems are more correlated than those with low CTv 
(< CGG/^gg > =2.14±0.37and < ^gg/^ss > =1.29±0.17, 
respectively). This evidence is in agreement with that 
found for clusters, and suggests a continuum of cluster- 
ing properties for all galaxy systems. 
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s (h"' Mpc) 

Fig. 6. The group-group (filled circles) and galaxy-galaxy 
(open circles) CFs, as computed for the volume-limited 
sample. 




s (h"' Mpc) 

Fig. 7. The ratio Ccc/Cgg^ the combined sample 
(thin line), and its nearby (dashed line) and distant (thick 
line) parts, as computed with the unweighted CF. 



In this context, it is appropriate to compare the groups 
of the volume-limited sample, characterized by the me- 
dian velocity dispersion — 214 km s~^ and by the 
mean group separation d ~ 16 h~^ Mpc, to rich clus- 
ters (cTv ~ 700 km s-i ; d - 50 h^^ Mpc; e.g. Zabludoff 



et al. |l993t Peacock & West |1992D . We fit ^gg(s) to the 
form ^(s) — {s/so)~'^ with a non-linear weighted least 
squares method and find 7 = 1.9 ± 0.7 and Sq = 8 ± 1 
h~^ Mpc. Note that groups show similar slope but signif- 
icantly smaller correlation length than optically or X-ray 
selected clust ers, for whic h sp > 15 h^^ Mp c (e.g . Bahcall 
& We st |1992| ; Croft et alJTgM^Abadi et al. |1998| ; Borgani 



et al. |199S| ; Miller et al. |1999D . Our resuhs agree with the 



predictions of those N-body cosmological simulations that 
also correctly predict the observed cluster-cluster CF (e.g. 
cf. our (so, d) with Fig. 8 of Governato et al. 1999). 



4-3. The unweighted CF 

In order to verify the stability of our results against varia- 
tions of the weighting scheme, we compute the unweighted 
CF, for the magnitude-limited sample. We find 

that, as in the weighted case, the amplitude of Cgg'('5) 
still significantly higher than the amplitude of S.^^ (s) , 
< Cgg /C^g^ > =1-18 ± 0.05. We also find that the ampli- 
tudes of ^qg'('^) ^^'-^ ^gg^ i^) ^^'^ both significantly lower 
than the weighted estimates. 

The differences between the results of the two weight- 
ing schemes rise from the fact that the weighted CF 
weights each volume of space equally and therefore bet- 
ter traces the clustering of more distant objects. In fact, 
when we divide the group/galaxy catalogs in two subsam- 
ples of equal size according to group/galaxy distances, 
the distant samples {V > 6680 km ) give CFs 

with higher amplitude, i.e. < CGG^distant/CcG^noarby > = 

2.15 ± 0.31 and < &ta„t/?g"g^earby >= 1-43 ± 0.08. 
Moreover, the distant samples {V > 6680 km s ^ ) give 



< e 



UW /tUW 
GG, distant/ '^gg, distant 



>= 1.43 ±0.12 in closer agree- 



ment with the result of the weighted analysis. The ratio 
^GG /^g^ for the whole sample, as well as for its nearby 
and distant parts, is shown in Fig. |^. 

As for a physical explanation, the fact that 
^gg^('^) <^ggi^) could be the consequence of a depen- 
dency of clustering on luminosity, since the unweighted 
CF estimator is more sensitive to the cl usteri ng of nearer, 
fainter groups/galaxies (e.g.. Park et al. 1994). In fact, the 
dependency of clustering on luminosity has been pointed 
out for the galaxy -galaxy CF (e.g_ 
Cappi et al 



1998 



19981) 



Benoist et al 
In addition 



1996 



the 



Willmer et al. 

greater strength of ^gg(s) could be explained by different 
clustering properties in dif ferent volumes of the Universe: 
e.g., Ramella et al. (1992) find that the strength of the 
galaxy CF is very high in the Great Wall. In the volume 
we examine the two biggest structures, the Great Wall 
and the So uthern Wall, both lie in distant regions (e.g. da 
Costa et al. 1994 ) and therefore their weight is larger in the 
weighted CF scheme. It is reasonable to expect also that 
distant groups, which are brighter (and presumably more 
massive) and which preferably lie in the two big structures, 
are more strongly correlated than nearby groups leading 
to the observed ^g^(s) <^gg(s) ■ 

5. Summary and conclusions 

We measure the two-point redshift-space correlation func- 
tion of loose groups, ^gg(s) , for the combined CfA2 
and SSRS2 surveys. Our combined group catalog con- 
stitutes the largest homogeneous sample available (885 
groups). We compare ^gg('S) with the correlation func- 
tions of galaxies, Cggis) , in the same volumes. 
Our main results are the following: 

1. Using the whole sample we find that groups are sig- 
nificantly more clustered than galaxies, < ^GG/^gg > 
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=1.64 ± 0.16, thus consistent with the result by 
Trasarti-Battistoni et al. (1997), based on a much 
smaller sample. This ratio can be considered a lower 
limit considering the possible presence of unphysical 
groups. 

Groups are significantly less clustered than clusters. 
In particular, we find 7 — 1.9 ± 0.7 and sq = 8 ± 1 
for 139 groups identified in a volume-limited sample 
{y < 7800 km s^^ , median velocity dispersion ^ 
200 km , and mean group separation d ~ 16 h~^ 
Mpc) . This result can be compared with that of galaxy 
clusters (sq ~ 15-20 h~^ Mpc for systems with ^ 
700 km s"i and d 50 h"^ Mpc; e.g., Bahcall & West 



1992; Croft et al. 1997; Abadi et al. 1998; Borgani et 



al. |1999| ; Miller et al. |1999|). 

3. There is a tendency of clustering amplitude to increase 
with group velocity dispersion Cv, which is the better 
indicator of group mass at our disposal. 

We conclude that there is a continuum of clustering prop- 
erties of galaxy systems, from poor groups to very rich 
clusters, with correlation length increasing with increas- 
ing mass of the system. 
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